Effect of annealing on the hyperfine interaction in InAs/GaAs quantum dots 
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The hyperfine interaction of an electron with nuclei in the annealed self- assembled InAs/GaAs 
quantum dots is theoretically analyzed. For this purpose, the anneaUng process, and energy struc- 
ture of the quantum dots are numerically modeled. The modeling is verified by comparing the 
calculated optical transitions with the experimental data on photoluminescence for a set of the an- 
nealed quantum dots. The localization volume of the electron in the ground state and the partial 
contributions of In, Ga, and As nuclei to the hyperfine interaction are calculated as functions of the 
annealing temperature. It is established that the contribution of indium nuclei into the hyperfine 
interaction predominates up to high annealing temperatures (Ta — 980 °C) when the In content in 
the quantum dots does not exceed 25%. The effect of the nuclear spin fluctuations on the electron 
spin polarization is numerically modeled. Effective field of the fluctuations is found to be in good 
agreement with the available experimental data. 

PACS numbers: 72.25.Rb, 78.67.Hc, 73. 21. La 



I. INTRODUCTION 

Hyperfine interaction of an electron with a nuclear spin 
ensemble is known to give rise to the most effective mech- 
anism of electron spin relaxation in quantum dotsis^i^ 
Due to limited number of nuclei in a quantum dot (QD) 
which interact with the electron spin, combination of the 
randomly oriented nuclear spins leads to non-zero total 
spin of the nuclear systemi^i^ This total spin acts on the 
electron spin as an effective magnetic field, B^^ with 
random magnitude and orientation. The electron spin 
rapidly precesses about this fluctuating field that results 
in a decay of the electron spin polarization in the QD 
ensemble. Typical times of the spin decay are of about 
a fraction of nanosecond for InAs/GaAs QDs.^ At the 
same time, the electron spin relaxation due to other pro- 
cesses such as electron-phonon interaction is a few orders 
of magnitude longer ^^li 

The hyperfine interaction strength depends on the 
number of nuclei covered by the electron wave function, 
thus on the electron localization volume.^ Scaling the QD 
effective size we can efficiently control the hyperfine in- 
teraction. A good way to change the QD size is the 
postgrowth thermal annealing the heterostructure with 
the QDs at relatively high temperatures. The annealing 
causes diffusion of indium from the QDs into the barrier 
layers that results in decreasing the potential well depths 
and its enlarging in real space. Besides, the annealing 
allows one to reduce the number of point defects and, 
thus, to improve the structure quality^ that suppresses 
the defect-related electron spin relaxation mechanism. 

In this paper, we theoretically consider the effect 
of postgrowth annealing on the electron-nuclei hyper- 
fine interaction in the self-assembled InAs/GaAs QDs. 
We offer a numerical model of the QD which is based 
on the real data obtained in experiments with such 



heterostructures.^'^^^^^ This model allows us to simulate 
the process of annealing and to calculate the distribu- 
tion profiles of In and Ga atoms over the heterostructure. 
Using this model we calculated the energy states of the 
carriers and the optical transitions in the annealed QDs. 
We fitted the parameters of our model using the experi- 
mental data on photoluminescence (PL) for a set of the 
heterostructures with the InAs/GaAs QDs annealed at 
different temperatures. This fitting allowed us to model 
the hyperfine interaction of an electron with the nuclei as 
a function of annealing temperature of the heterostruc- 
ture. Finally, the result of these calculations made it 
possible to quantitatively describe the electron spin de- 
polarization via the hyperfine interaction with nuclei and 
suppression this effect by a longitudinal magnetic field. 



II. QUANTUM DOT MODEL 

Typically the self- assembled InAs/GaAs quantum dots 
have the truncated pyramid or lens shape with base diam- 
eter d= 15-30 nm and height h = 5-15 nmJ^ The size 
of the QDs depends on the growth parameters, in par- 
ticular on the nominal thickness of the deposited layer of 
indium. Besides, there is a spread of sizes of QDs in the 
QD ensemble. To be specific, we consider heterostructure 
investigated in Refs. [sI ITqIITi ] . Its cross-section image ob- 
tained by scanning electron microscopy (SEM) is shown 
in Fig. [TJ Though the spatial resolution is not high we 
can estimate the base diameter of the QDs to be of about 
20-25 nm. A higher spatial resolution can be obtained 
by transmission electron microscopy, but such data are 
not available for this heterostructure. 

Therefore we use indirect data on the QD shape and 
size which can be extracted from the analysis of the PL 
spectra. As it is well known, the quantization energy 
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where x{r,t) is a position-dependent function of indium 
fraction in the In^^Gai-a^As sohd solution forming the 
QD annealed during the time and D is the diffusion 
constant. The diffusion equation ([T|) should be supple- 
mented with the initial condition on the function x. In 
our model we assume that the unannealed QD is pure 
InAs, so that the initial conditions are: 



FIG. 1: SEM image of the cross-section of the heterostructure 
with unannealed InAs/GaAs QDs. 



of carriers and, therefore, the energies of optical tran- 
sitions depends on size and shape of a QD. We choose 
these parameters to obtain correspondence between the 
experimentally observed and calculated energies of op- 
tical transitions between the lowest and excited states. 
The lowest optical transition depends mainly on the size 
of the QD. However, the energy distance between the 
lowest and excited states and relevant optical transition 
energies depend on the ratio of the height to base diame- 
ter of the QD. We use these circumstances to determine 
separately the height and the base diameter of the QD 
under study. We ignore the statistical spread of sizes of 
the QDs in ensemble. Moreover, for simplification of cal- 
culations we consider the QD with cylindrical symmetry 
and smooth bell-like shape. The QD model is schemati- 
cally shown in Fig.[2l The height of the QD Hqd = S nm 
and the base diameter (at height OAhgo) dgo = ^0 nm. 
We put the QD on a thin layer of InAs {hwL = 0.283 nm) 
to model the wetting layer (WL) which is inevitably ap- 
pear when growing the QDs in the Stransky-Krastanow 
growth mode. 



III. GALLIUM AND INDIUM 
INTERDIFFUSION PROCESS DUE TO 
ANNEALING OF THE HETEROSTRUCTURE 

Postgrowth annealing of a heterostructure with self- 
assembled InAs/GaAs quantum dots leads to indium and 
gallium inter diffusion.-j^ii^ii^ Like to other authors, we 
consider this process in a model of continuum because the 
effective scale which we are interested in is much larger 
then the lattice constant. Besides we assume the diffusion 
coefficient to be independent of space coordinates. So we 
describe the diffusion by the Pick's law: 



9x(r, t) 
dt 



x(r,0) 



1, inside QD/WL 

0, inside the barrier layers. 



(2) 



In the accepted model of cylindrical symmetry of the 
QDs, we use the cylindrical co-ordinates. Defining axial 
axis z to be perpendicular to the plane of the quantum- 
well layer (WL), and radial axis p that lies in the 
quantum- well plain (see Fig. [2]). In these co-ordinates, 
the function x(r^t) has got a separable form in azimuthal 
angle x{r,t) = xi^i P^t)^{(p). Because the cylindrical 
symmetry is conserved during the annealing process, the 
indium distribution x does not depend on azimuthal an- 
gle, i. e. ^{(p) = 1. Dividing Eq. ^ on D and separating 
the variables, we can write the equation for the function 
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(3) 



The initial conditions for the function x coincide with 
those for x [see Eq. ([2])]. 

We choose the computation region to be a cylinder 
with the height, Hcyi = 80 nm, and the diameter, 
Dcyi = 80 nm. As we found, these values are much larger 
then the diffusion length of indium atoms in InAs/GaAs 
semiconductors, < 10 nm, at annealing temperatures up 
to 1000 ^C. For the rigorous description of the prob- 
lem, Eq. (|3]) should be supplemented with appropriate 
boundary conditions. The Dirichlet boundary condi- 
tions, X = 0? ai'c imposed on boundaries 1 and 3 because 




GaAs 
barriers 



DAx{r,t) = 0, 



(1) 



FIG. 2: Model of the quantum dot: cross-section of the het- 
erostructure and the computational region (see text for de- 
tails) . 
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these boundaries are far from the QD (see Fig. [2]). We 
also assume that the indium flux through boundary 2 is 
zero because, in the neighborhood of this boundary, the 
In atoms diffuse from the WL in the perpendicular di- 
rection which is evident from the local symmetry of the 
problem. This assumption can be described by the Neu- 
mann boundary condition, n ■ Vx = 0, where n is the 
outward normal vector to the boundary. On boundary 4, 
the Neumann boundary conditions also must be imposed 
for nulling the diffusion flux across the symmetry axis of 
the dot. 

Using finite element technique, we solve the described 
diffusion problem with different diffusion constants and 
the fixed annealing time interval, t = 30 sec. We assumed 
the Arrhenius equation for the temperature dependence 
of the diffusion constant: 



D{Ta)=Doexi> 



(4) 



where is the annealing temperature, Ea is the activa- 
tion energy of the interdiffusion process, k is the Boltz- 
mann constant, and I^o is the pre-exponential factor. We 
used Ea and I^o as fitting parameters whose values were 
optimized to get the best agreement between the experi- 
mentally measured and calculated PL spectra of the an- 
nealed QDs (see Sec. lTVC]) . We found that Ea = 1.23 eV 
and 1^0 = 8.5 x 10~^^ m^/s. The obtained value of Ea 
is approximately three times smaller than reported one 
in Ref. [13] for annealed InGaAs QDs. The diffusion 
length Ld = y/[Dt] , which is usually measured in exper- 
iments, is also larger in our calculations {Ljj = 3.6 nm 
vs Ld ^ 1.5 nm in Ref. [9] at 900 °C). The origin of 
discrepancy of these quantities is unclear. We found that 
variations of parameters of our QD model (size of QD, 
band-offset, strain energy, see Sec. llVCj) in the ranges 
when the calculated optical transitions agree to experi- 
ment (see Sec. lIVCp do not change the activation energy 
noticeably. 

Fig. [3] shows examples of the In distribution over the 
heterostructure calculated for the QDs unannealed (a) 
and annealed at different temperatures (b-d). As seen, 
annealing of the heterostructure leads to the rapid disso- 
lution of the InAs QD into the GaAs barrier layers. The 
average indium content does not exceed 25% for QDs an- 
nealed at temperatures above 980 ^C. Besides, the QD 
volume increases with annealing temperature. 



IV. ENERGY STRUCTURE OF THE 
ANNEALED QUANTUM DOTS 

A. Potential profiles and effective masses 

Diffusion of indium from the QDs into the barrier lay- 
ers leads to modification of the profiles of the valence 
and conduction bands. Using the calculated indium con- 
tent x{r) for the annealed QD, we modeled the three- 



dimensional potential profiles as well as the profiles for 
the effective masses of carriers. 

For this purpose, we used a linear approximation for 
the bandgap of In^^Gai-ajAs using relevant quantities for 
InAs and GaAs:^^ 



;iInAs 



x{r) + E\ 



GaAs 



[l-x(r)], (5) 



where E^^^^ = 0.415 eV and E^""^^ = 1.519 eV are the 
bandgaps for InAs and GaAs, respect ivelly. Besides, we 
used the band-offset ratio, Qe/Qh = 76/34, along; with 
the the data for InAs/GaAs taken from Refs. [l6lfl7l|. 
There is, however, an important problem which compli- 
cates this point. Namely, the large mismatch between 
the InAs and GaAs lattice constants gives rise to the 
large built-in strain which considerably affect the poten- 
tial profiles^aiO^aa^^^ 

We include this strain in our model in a simple way. 
At the begining, let us consider the unannealed QD. The 
average value of hydrostatic strain in the InAs/GaAs 
QDs, 7-10 %, which obtained in calculations by differ- 
ent authors using continuum elasticity theoTjr^^^^ 
and valence-force- field Keating model, ^^^—^^i^ leads to 
decreasing the well depth in the conduction band by 350 - 
500 meV. Moreover, the strain weakly changes over the 
QD volume and, therefore, the constant potential ap- 
proximation gives rise to the good results. We assume 
in our model the strain-induced energy shift of the con- 
duction band by = 450 meV, so that the electron 
confining potential is 390 meV. Using this value, we are 
able to correctly describe the optical transition energies 
observed experimentally (see Sec. IIVC]) . 

The confinement potential shape in the valence 
band is more complex than that in the conduction 
band.^'^^5^i2£i^'^^ However, Califano et al?^ used a simple 
model of the constant potential for the valence band and 
determined the potential well depth from comparison of 
their calculations with the experimentally found energies 
of optical transitions. Using similar method, we obtained 
the strain energy in the valence band = 90 meV which 
corresponds to the hole confining potential 175 meV. 

The annealing leads to redistribution the indium con- 
centration and, as we believe, to relaxation the built-in 
strain. We assume the strain energy linearly depends on 
the indium concentration in the annealed QDs. Thus, 
using the calculated distribution x(r), the bandgap (see 
Eq. ([5])), the band-offset ratio, and the linear strain de- 
pendence on x(r), we can determine the potential profiles 
by means of equation: 

Ve,h{r) = Qe,h (^g(GaAs) - Eg{r)) - ■ x{r). (6) 

The cross-sections of potential profiles along z-axis thus 
determined are shown in Fig. [H 

The strain also affects the effective masses of carriers. 
In the dot material the compressive stress alters the cur- 
vature of the bulk bands, causing the effective masses 
to differ from the unstrained ones. We have used the 
values of electron effective masses in the GaAs barriers 
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FIG. 3: (Color online) (a-d) Cross-section of calculated distributions of indium content x over the heterostructure unannealed 
and annealed at different temperatures. Highest values of indium concentration are marked in percents for each QDs. Thin 
black line shows the unannealed QD shape, (ae-de) Distributions of electron density for the ground state in the unannealed 
and annealed QDs. The cross-sections of probability density isosurface at level O.33|'0|max are shown by thin yellow ellipse, 
(ah-dh) Distributions of the ground state hole densities. The notations are similar to those used for the electron density 
distributions. 



and in the strained InAs QD, m*(GaAs) = 0.0665mo 
(Ref. [H), and m*(InAs) = 0.04mo (Ref. [23]), respec- 
tively. The hole effective masses used in our calcula- 
tions are: m^(GaAs) = 0.3774mo, and mJ^(InAs) = 
0.341mo.^^ Here mo is the electron mass. In the an- 
nealed QDs, we used the linear dependence of the carrier 
effective masses on the indium concentration: 



^ir) 



, (InAs) •x(r ) +m: . (GaAs) • [1 -x(r)] . (7) 



B. The electron and hole energy states 

To compute the electron and hole states for 
the annealed heterostructure, we solve the one-band 
Schrodinger equation in the effective mass approxima- 
tion: 



_m*(r) 



V?A(r) 



V{r)tP{r) =Ei){r), (8) 



where fi is the Planck's constant divided by 27r, m*{r) 
and V{r) are position-dependent the potential energy 



and the electron (or hole) effective mass [see Eqs. (j6l 
E is the carrier energy, and '^(r) is the envelope wave- 
function. Because the annealed QD/WL system in our 
model has a perfect cylindrical symmetry, we again can 
write the carrier position in the cylindrical co-ordinate 
system, r = (z^p^cp). We can partially separate the co- 
ordinates in the total wave-function, ^^(r) = 77(2;, p)Q{(p), 
and rewrite the Schrodinger equation in the cylindrical 
co-ordinates as: 



2 



d_ 

dz 



1 df] 

m* dz 



p dp \m* dp 

'. ^ 

2m* dcp'^ 



e- 



VrjS = ErjS. 



(9) 



Dividing this equation by r]{z^ p)Q{(p)/ {m^ p'^) and rear- 
ranging its terms, we come to the two independent equa- 
tions: 
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dz 



1 drjn 
m* dz 
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p dp \m* dp 



(10) 
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2 e 9(^2 2 

The last equation for Q{(p) can be solved analytically: 

1 



(11) 



exp [iriLp] . 



(12) 



where n should be an integer to get one- valued function 
[e(27r) = 9(0)]. 

Then, we should discuss the physically substantiated 
boundary conditions for function 7^^(z,p). We are only 
interested in electronic (hole) states confined in the QD. 
Therefore, we can impose the Dirichlet boundary condi- 
tions on boundaries 1-3, r]n{z^p) = 0, to provide the 
wave- function dumping in the barrier layers (see Fig. [2]). 
Then we consider the boundary conditions on boundary 
4 (symmetry axis). There are two different types of the 
conditions. When n ^ 1, the Dirichlet boundary con- 
ditions must be imposed to ensure that the last term 
in Eq. (p!Q|) does not diverge on symmetry axis, p = 0. 
When n = 0, we employ Neumann boundary conditions, 
^ • ^Vn{z^p) = 0, to ensure VV^ existing. Here n is 
the outward normal vector to the boundary. Apart from 
these assumptions, we must impose constraints on func- 
tion r]n{z,p), 
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\rin{z,p)\'^pdzdp = 1, 



(13) 



to satisfy the normalization conditions for the wave- 
function. 

Using finite element technique, we solve eigenvalue 
problem described by Eq. (p!Q|) for electrons and holes in 
the QDs, unannealed and annealed at different tempera- 
tures. The annealing leads to a modification the carrier 
density shape and to an increase the localization area 
at high annealing temperatures. Namely, at low anneal- 
ing temperatures (less than 820 ^C), the shape of the 
QD is changed a little that almost does not affect the 
carrier clouds. However, the changes of potential well 
shapes promotes some additional localization of the car- 
riers. The annealing at greater temperatures leads to 
spherical-like shapes of the carrier density distributions, 
which corresponds to the spherical-like shape of the in- 
dium concentration distribution. These effects are illus- 
trated in Fig. [31 where the electron densities (ae-de) and 
hole densities (ah-dh) for the lowest states in the unan- 
nealed and annealed QDs are shown. 

Also the annealing of the QD leads to an increase in the 
electron and hole energy levels. The latter effect is illus- 
trated in Fig. m As seen, the potential well depth is de- 
creased that provides the decreasing the energy distance 
between the energy levels and bottom of the potential 
well. However, the energy gap between the electron and 
hole states increases. This affects the energies of optical 
transitions in the annealed QDs. 
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FIG. 4: (Color online) Potential profiles for the valence and 
conduction bands along the symmetry axis of problem (/> = 0) 
for the QDs unannealed and annealed at different tempera- 
tures. The energy position of the lowest electron and hole 
states are shown by thin horizontal lines. The calculated en- 
ergies of the ground optical transitions are shown in eV. 



C. Optical transitions and comparison with PL 
spectra 

The next step of modeling is describing the optical 
transitions in annealed QDs. To calculate their energies, 
we should take into account the energy of the electron- 
hole Coulomb interaction. Rigorous solution of this prob- 
lem is extremely difficult because we must consider the 
quantum mechanical problem for electron-hole pair in the 
configuration space with six degrees of freedom. How- 
ever, since the localizing potential for carriers in a QD 
is much larger then the potential due to the Coulomb 
interaction, we may ignore the Coulomb correlations in 
the motion of the electron and hole in the QD, and cal- 
culate the Coulomb interaction energy, Ec , for two fixed 
spatially distributed charged clouds:— 



km 



II 



We - rh\ 



d^red^Vh, (14) 



where e is the elementary electronic charge, Sq is the vac- 
uum dielectric constant, e is the the average dielectric 
constant in the annealed QD, rg, and Vh are the electron 
and hole positions, and (rg), and iljj^{rh) are the en- 
velope wave-functions of electron in the kth energy state 
and hole in the mth energy state, respectively, which are 
calculated above in the framework of the one-band prob- 
lem. Calculations show that is of about 21 meV for 
unannealed QD and decreases down to 15 meV for the 
QD annealed at 980 ^C. This energy is slightly smaller 
when the carriers are in the excited states. 

Then the energies of optical transitions are calculated 
from the simple equation: 



Ekr. 



E^ -^EY'^E^-E, 



km 
C 5 



(15) 



where E^ and E]^ are the electron and hole energies with 
respect to lUa^Gai-a^As conduction band bottom and va- 
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FIG. 5: Color online) Calculated energies of the lowest optical 
transitions (solid vertical lines) in comparison with the PL 
spectra for the QD's unannealed and annealed at different 
temperatures. The dashed line shows the energy position of 
the optical transition between the first excited electron and 
hole states in the unannealed QD. 



le l/i, 2e 2/i, etc. This interpretation also con- 
firmed by the available experimental data obtained at 
the high excitation power in magnetic field up to 28 T 
which are not show here. 

As seen from the Fig. [5l the calculated transitions 
Oe Oh well reproduce the experiment for all the series. 
It should be noted, however, that the very good agree- 
ment in energies of the calculated and measured tran- 
sition le Ih is obtained for the unannealed sample 
only. For the annealed QDs, the calculated energy of the 
transition is larger by about 20 meV which is compara- 
ble with the distance between the Oe Oh and le ^ Ih 
transitions. This disagreement probably related to over- 
simplifying the strain distribution in our model. How- 
ever, because we only interested in the lowest electron 
states which are involved in the spin relaxation problem 
under consideration, this disagreement appears to be not 
very important. 



V. HYPERFINE INTERACTION OF THE 
LOCALIZED ELECTRON WITH NUCLEI 

A. Effective magnetic field of the nuclear spin 
fluctuations 



lence band top, respectively. Eg is the lUa^Gai-a^As band- 
gap energy at the QD center. 

We calculated the optical transitions between the low- 
est electron and hole energy states as well as between 
the excited states. We found that the energy of the 
lowest optical transition mainly depends on size of the 
QD. At the same time, the distance between the lowest 
and excited transitions is mainly determined by the ra- 
tio height /diameter and the shape of the QD. The strain 
energies for valence and conduction bands, which we also 
used as fitting parameters, mainly influence on the shift 
of optical transition series, and weakly affect the distance 
between the transition energies. These facts allow us 
to uniquely determine the height and the diameter of 
the QD by comparing the calculated and experimentally 
measured PL spectra for the unannealed sample. Then, 
the annealing was modeled as it is described below in 
Sec. mil We used the interdiffusion activation energy 
Ea = 1.23 eV and the prefactor Dq = 8.5 x lO"^"^ m^s 
[see Eq. @] as fitting parameters for all the series of 
annealed samples. 

For comparing with the calculations we measured the 
PL spectra of a set of the samples using a standard ex- 
perimental technique. The spectra are measured at ex- 
citation photon energy E^x = 2.54 eV at temperature 
T = 1.6 K. They are shown in Fig. [5l It is experimen- 
tally and theoretically founded that, the first PL peak 
corresponds to recombination of the electron and hole in 
the lowest states, Oe Oh. The next PL peaks corre- 
spond to recombination of an electron and a hole which 
are in excited states with the same quantum numbers: 



As it is discussed in Introduction, the electron spin po- 
larization is efficiently destroyed in QD ensemble via hy- 
perfine interaction with randomly oriented nuclear spins. 
Theoretical justification of the electron spin relaxation 
mechanism in QDs was reported in Ref. [^]. The general 
idea of this mechanism can be described quantitatively 
as follows. The interaction of the electron and nuclear 
spins is determined by their hyperfine Fermi interaction: 



(16) 



where S, Ij are electron and jth nucleus spins; = 
[WiTjiiBl-i^j /i^IjVo)] • \uc{Rj)\'^ is the hyperfine coupling 
constant with the jth nucleus; /i^ is the Bohr magne- 
ton; ^0 is the unit cell volume; /j, /ij and Rj are spin 
magnitude, magnetic moment, and position of the jth nu- 
cleus; il^{Rj) and Uc{Rj) are the electron envelope wave- 
function and the Bloch function at the nuclear site. 

Due to the limited number of nuclei in a QD which in- 
teract with the electron spin, random orientation of the 
nuclear spins give non-zero total spin which has a magni- 
tude fluctuating from dot to dot. The total nuclear spin 
acts on the electron spin as an effective nuclear hyperfine 
magnetic field, Bn. We consider the non-polarized and 
not interacting nuclear spins, since the magnitude and 
orientation of the effective field are random and can be 
described by the normal distribution^^ 



wb 



{V2nABy 



■ exp 



2AI 



(17) 
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with variance determined by: 



(18) 



where bj = [vq / {jiiBge)]^^ is the effective mag- 

netic field of a single nuclear spin acting on the elec- 
tron. Here ge is the electron g-factor. Unlike paper by 
Merkulov et al.^^ we defined variance [see Eq. (p!8|) ] so 
that it approximately corresponds to half width at half 
maximum of normal distribution given by Eq. (pT|) . 

In the InGaAs QDs, there are three types of nuclei. 
We consider the fluctuating field as a sum of three inde- 
pendent contributions with normal distribution of each 
of them. The total variance squared known to be calcu- 
lated as a sum of variance squares of variate independent 
contributions, and, therefore, we can write: 



Ga\2 



(A^^)2, 



(19) 



where A^, A^^, and A^^ are the partial contributions 
of each types of nuclei: 



l)h]^, ^ = In, Ga, As. (20) 
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Here the sum goes only over each types of nuclei in crys- 
tal lattice. Supposing that the electron envelope wave- 
function is constant over the crystal unit cell and replac- 
ing the sum over unit cells by the heterostructure volume 
integration, we obtain: 



{^'Sf = Ci.j\^{rtx{r)d\, 
{A%n^=CGaJmrt[l-x{r)]d' 



(21) 



where x{r) is the indium fraction in Ina;Gai_a;As solid so- 
lution forming the QD, and constants are determined 
by: 

1 

C^ = -k{h^^)7 ^2^0, e = In, Ga, As. (22) 

Actually x{r) determines the probability to find the in- 
dium nucleus in the position r. 

Next, we introduce the new parameter, the average 
effective indium fraction in the QD: 



j\ilj{r)\^x{r)d^r 



(23) 



j\tP{r)\HH ■ 
Using the electron localization volume, defined in Ref. [4] : 



Vl 



\^{r)\'d^r 



(24) 
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FIG. 6: (Color online) (a) Dependencies of the electron lo- 
calization volume (red) and of the average effective indium 
fraction in the QD (black) on the annealing temperature; (b) 
The total nuclear field variance and the partial contributions 
of each type of nuclei as a functions of the QD annealing 
temperature. 



we come to the final expression for the total variance of 
the effective nuclear field: 

^l = i- (^InX + CGa(l - X) + Cas) • (25) 
VL 

We calculated the effective indium fraction and the 
electron localization volume as a function of the QD an- 
nealing temperature. Then, using Eq. ([25|) . we calculated 
the total variance and the partial contributions of each 
types of nuclei. In these calculations, we used the semi- 
conductor parameters taken from Ref. [26^ : 



Nuclei 


In 


Ga* 


As 


Nuclear spin / 


9/2 


3/2 


3/2 


Hyperfine constant A (/ieV) 


56 


42 


46 



* Average between ^^Ga and '^-'^Ga. 

Results of the calculations are shown in Fig.[6l As seen, 
the electron localization volume increases with tempera- 
ture which is due to decrease of the localizing potential 
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depth in QD and the increase of the QD size. Also, the 
dissolution of the QD in the barriers leads to decreasing 
of the effective indium concentration x. Two these ef- 
fects result in the decrease of the effective nuclear field 
variance. From physical point of view, the variance de- 
crease is due to (i) decrease of contribution of each nu- 
cleus to hyperfine interaction because of reducing of the 
electron density at the nucleus and (ii) averaging of the 
contributions over increasing number of nuclei. As seen 
from Fig. [6] (b), the gallium contribution to the disper- 
sion slightly increases with annealing up to temperature 
Ta = 850 which is caused by the gallium diffusion 
into the QD. However, this effect does not influence upon 
behavior of the total variance because the indium contri- 
bution is dominating at all annealing temperatures due 
to the large indium nuclear spin /in = 9/2. 




W 0.0 



-200 -100 100 200 

Magnetic Field, B (mT) 



B. Suppression of the nuclear spin fluctuations by 
the external magnetic field 

The electron spin relaxation caused by the nuclear fluc- 
tuating field can be suppressed by applying the external 
magnetic field. ^ In the presence of the external magnetic 
field, the electron spin precesses about the total field, 
= ^ext Bjsf (see inset in Fig. [7|). At sufficiently 
large external field, the nuclear spin fluctuations almost 
do not contribute to the total field and, therefore, the 
electron spin polarization does not decay. 

Let us analyze this effect in more detail. We con- 
sider behavior of the projection of the electron spin {z- 
projection) onto the axis of the optical excitation which 
is typically measured in experiment. The direction of the 
external magnetic field (longitudinal field) coincides with 
z-axis. The electron spin z-projection should be averaged 
over many periods of the spin precession about the to- 
tal field and over the QD ensemble. The time averaging 
allows us to calculate the constant component measured 
in such experiments. Although all the QDs in the en- 
semble are identical in our model, the averaging over the 
ensemble arises because of the random magnitude and 
orientation of the nuclear spin fluctuations. The spin 
projection is calculated as: 



< >-- 



S,{B^^t)wB{BN)d?B, 



N- 



(26) 



Here, the probability distribution wb{Bm) is given by 
Eq. (HZD and ^^(Bext) is: 



5,(Bext) = 5o[cOS^(Bext)]' = 

{Bnz + -Bext)^ 



So 



(27) 



Nn 



where Sq is the initial electron spin polarization (we as- 
sume ^0 = 1/2), and is the angle between z-axis and 
the total magnetic field direction (see inset in Fig. [7j)- 
Note, that Eq. ([26|) is the averaging over the ensemble 



FIG. 7: (Color online) Calculated external magnetic field de- 
pendencies of electron spin polarization in the QDs unan- 
nealed and annealed at different temperatures (solid lines) 
comparing with experiment^^ (circles). Dashed line is the ap- 
proximation of calculations by the Lorentz function [Eq. (|3Qp ] 
for the unannealed sample. Right inset shows the mechanism 
of suppression of the nuclear spin fluctuation effect by an ex- 
ternal magnetic field. 



of the electron spin projection and Eq. (|27|) is a result of 
the time averaging of the electron spin precession. Fast 
precession of the electron spin about 5t conserves only 
the projection of the initial spin 5*0 onto the direction 
of B^ so that S = So cos 0. Therefore, the measurable 
quantity is 5'^ = 5* cos 6 = So cos^ 6. 

Using the values of found above [see Fig. [6] (b)], 
we calculated the dependencies of Sz as functions of the 
external magnetic field for the QDs annealed at differ- 
ent temperatures. The results are shown in Fig. [71 (solid 
lines). As seen, each calculated dependence has a dip 
around zero external field which is due to the depolar- 
ization of the electron spin by the effective nuclear field. 
This depolarization is incomplete, namely (l/3)5'o is still 
conserved. The result can be understood from a simple 
consideration. We may replace the arbitrary oriented nu- 
clear fields by three components oriented along the 7/, 
and z co-ordinate axes with equal probabilities, P = 1/3, 
and write the magnetic-field-dependent value of the elec- 
tron spin upon action of these components: 



Bn \\x,y: 
BnW z: 



Sz = So 
Sz = So- 



^ext 



BL^^Bl 



N 



(28) 



As seen, the nuclear fields along x and y axes totally 
depolarize the electron spin at zero external magnetic 
field, and the nuclear field along z-axis keeps its initial 
value. Assuming that all the components of nuclear field 
have identical average variances, which can be defined as 
< Bj^^ y z ^~ (2^^)^' average electron spin polar- 
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ization can be written: 

Rearranging terms in this equation gives the magnetic 
field dependence with the Lorentz-like dip: 

'•^"°('- 1 + [B.!/(2Ab)p )^ <'°' 

As seen from Fig. [71 the calculated field dependencies of 
< 5'^ > can be perfectly approximated, in the average, 
by this Lorentz function. 

The dip widths at half minima are of several tens of 
milli-Tesla and decrease with the annealing temperature. 
That refiects the decreasing of the hyperfine interaction 
strength in conformity with the above discussion. We 
can characterize the interaction strength by the averaged 
nuclear field, — 2Ab, which we define as a half width 
at half minimum (HWHM) of the dip. This means that 
and its dependence on the QD annealing temperature 
can be obtained from those quantities for shown in 
Fig. El (b). 

The electron spin polarization in the QDs annealed at 
900 was experimentally studied in the paper by Cher- 
bunin, et al.^ They have measured the circular PL po- 
larization of the singly negatively charged QDs and have 
found that the polarization is closely related to the spin 
orientation of the resident electrons. The magnetic field 
dependence of the PL polarization is found to reveal a 
dip around 5ext = 0. One of the experimental curves 
measured at the excitation density 4.4 W/cm^ is shown 
by circles in Fig. [71 As seen, the curve is very similar to 
those calculated theoretically. A small disagreement is 
partially related to the experimentally found dependence 
of the HWHM on the excitation density.^ With decreas- 
ing the excitation density down to P = 0.5 V/cm^, the 
HWHM increased up to approximately 30 mT. The cal- 
culated value of HWHM, Bn - 27 mT, for = 900 
is close to experimental one. 



VL CONCLUSION 

Theoretical modeling of InAs/GaAs QDs allowed us 
to simulate the effect of nuclear spin fiuctuations on the 
electron spin polarization. We determined the electron 
localization volume and the effective indium fraction in 
the QDs for different annealing temperatures. Due to 
interdiffusion of In and Ga during the annealing pro- 
cess, the QD size increases and, correspondingly, the 
electron localization volume considerably increases (from 
~ 1700 nm^ for unannealed QD up to ~ 4900 nm^ for 
the QD annealed at 980 ^C). At the same time, the dis- 
solution of the QD results in decreasing the effective in- 
dium concentration in the QD. We calculated the partial 
contributions of the indium, gallium, and arsenic nuclei 
to the effective magnetic field of the nuclear spin fiuc- 
tuations and found that the hyperfine interaction is de- 
termined mainly by the indium contribution. The effect 
of the fiuctuations decreases with the annealing tempera- 
ture due to (i) the increasing number of nuclei interacting 
with electron and (ii) the decreasing indium concentra- 
tion. The average magnitude of the effective hyperfine 
field decreases from ~ 41 mT down to ~ 15 mT with the 
annealing temperature up to 980 ^C. Finally, we modeled 
the suppression of the nuclear spin fiuctuation effect in 
the longitudinal magnetic field. The calculated dip of the 
electron spin polarization is very similar to that observed 
in the experiment.^^ 
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